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Excitons have recently been shown to deviate from pure bosons at densities a hundred times 
' smaller than the Mott density. The corresponding calculations relied on the unscreened excitonic 

I ground state wavefunction. A consistent inclusion of screening, by use of the fundamental eigen- 

. function of the Hulthen potential, vindicates this approximation. 

PACS numbers: 71.35Lk 

■ I. INTRODUCTION 

(N 

In a recent work [Q, the following question was addressed: Up to what density can excitons be treated as bosons? 
When the density of electron-hole (e-h) pairs in semiconductors is very large, the carriers form an e-h plasma, so 
^ ' that excitons do not exist anymore. From dimensional arguments, we could expect the cross-over to appear when the 
distance between two excitons is of the order of their size, i.e., when XNa^/V ~ 1, a-x being the exciton Bohr radius, 
^ N the number of excitons in the sample volume V, and A a dimensionless factor. The high-density regime in which the 
Coulomb energy is dominated by the kinetic energy is valid when the parameter r^, defined by TV 1 7r(rs a^)^ /V = 1, 
""1 is small compared to 1. This leads to A = 47r/3 ~ 4. The prefactor A can also be obtained by evaluating the carrier 

■ density for which, due to Coulomb screening, no excitonic bound state survives. The accepted result ^ for this Mott 

a density leads to A ~ 1. 
I ' Still, excitons are not perfect bosons because they "feel" each other not only through Coulomb interaction but also 
through Pauli exclusion between their electrons and holes [§,1). More precisely, the consequences of the fcrmionic 
^ character of excitons induced by Pauli exclusion between their constituents lead to prefactors A « 100 j^]. This final 
O result indicates that our calculations using the three-dimensional fundamental wavefunction of the unscreened exciton, 
which admittedly lack proper consistency, should not be too far from the truth, because screening effects are small in 
that density regime. 

In the present Letter, we justify this assumption by considering the fundamental eigenfunction of the Hulthen 
^ ' potential [||-^, which has been widely used to describe three-dimensional screening. The ensuing calculations are 
analytical throughout, and we basically obtain the same two orders of magnitude for A as previously. If anything, 
taking account of screening slightly increases A. 



o 



II. NOTATIONS 



Let us briefly summarize the key quantities introduced in Ref. For the sake of simplicity, we will restrict again 
to excitons with zero total momentum. In terms of the free electron and hole 6^^, creation operators, the exact 
^ [ creation operator for the ground-state exciton reads 

=J2 4>kalbl,, (1) 

O . k 
u ■ 

^ ] where (f>k is the ground-state exciton wavefunction in momentum space. From Eq. (P, we get 
•I— I . 

X; D = l-[B.B^=J2\^kmalak + bl,b-k)- (2) 

Obviously, D vanishes for perfect bosons. 

The relevant state in our problem is the A^-pair state constructed from the ground-state exciton 
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This state is the valid picture for the ground state of N pairs in the low-density limit, as shown by Keldysh and 
Kozlov We also showed that the quantity Fpf defined by 

{0\B^ {B^f\0) =NIFn (4) 

plays a central role since it appears in all expectation values taken on j^'^^^). If B were a perfect boson operator, we 
would get F;v = 1. Actually, the general expression of F^ reads |l|,p^ 

m=l ^ ' 

where the [im] are determined by ^ mim = N , after we have set 

a„, = Y,\hr- (6) 

k 

From the abovementioned expressions we could derive the main result of Ref. namely that excitons can be 
considered as true bosons provided that the expectation value of D defined in Eq. (H) is much smaller than unity, 
which translates into 

2 7Vcr2<l. (7) 



III. HULTHEN FUNDAMENTAL WAVEFUNCTIONS 

In contrast to Rcf. [Q , we now want to include some measure of screening in the wavefunction cjjk ■ This is possible 
by using the Hulthcn potential P-|9|Jl^ 

, , 4i? 1 

Vnir) , (8) 

g e~ - 1 

where R is the three-dimensional exciton binding energy, its Bohr radius, and g a dimensionless constant inversely 
proportional to the screening constant derived from the Debye-Hiickel or Thomas-Fermi models. Only a finite number 
of bound states remains; they all disappear for g < 1. Vh has an interesting property: all s eigenfunctions are known 
analytically. In the present study however, we are interested in the lowest-lying bound state only. In real space, its 
eigenfunction reads g 

r 

1 1 2 I T~ sinh(^— ) 

The analytical expression of (pk can thus be deduced by a mere Fourier transform of Eq. (|^): 



_ / 647ra3 1 1 

i V ^g^> [(1 - 1)2 [(1 + 1)2+^^2]' i^^i 

where K = ka^- In the g — > +oo limit, we recover the expected result — z/{l + if^)^, where z — GAira^/V. 



IV. CALCULATION OF THE NEW CRITERION 

By definition, the tJm's given by 



2 7r2a3 V y ^ h [(1 - i)^ + /f2]2m [(l + 1)2 + ;^2]2m 
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are now functions of the screening parameter g. The first ones are respectively 



ai(g) = l (12) 

, , 33 2 1 5 1 \ 

"^(^) = T28" 1^-3^+337; 

4199 ^2 A 12 1 6 1 28 1 63 1 \ 
0-3 (.9) - 32768 ^ 1^ ~ 13 52 + 13 ^ ~ 221 ^ ^ 4199 g8 j 

334305 .3 / 18 1 225 1 100 1 675 1 594 1 143 1 \ 
0-4(5) - 4194304^ V ~T7g2 + 323^~32356'^ 7249 ~ 37145 ^ + 111435 'g^ J 



with 



Z — — 

Actually, the general result can be expressed analytically as 



(13) 
(14) 
(15) 

(16) 



--^^^-('y^-' 1^1^ ^(2-2.n,|2.n+i;^), (17) 

where 2F1 is the hypergeometric function; when g —* +00, this function goes to 1, while z goes to z, so that we recover 
the result of Ref. ||l| for the unscreened exciton. 

Note that g appears at different places in Eq. (p^). The first effect of screening amounts to a "renormalization" of 
the excitonic Bohr radius from Qx to ax /{I — Because of the decreased electron- hole interaction, the exciton 

binding energy should decrease and the Bohr radius should therefore increase, tending to infinity at the dissociation 
limit. This is an unsurprising behavior. None the less, g also appears in the argument of the hypergeometric function, 
so that the simple physical picture of a dilated exciton is not strictly true. 

From Eq. (Q), we finally obtain for our perfect boson-exciton criterion 

1 - - — + — — 

2Na,ig) = 33n ^^l^'yj^l^' (18) 



V. DISCUSSION 



Let us now discuss how our consistent introduction of screening affects the criterion for a perfect boson behavior of 
excitons. We have A equal to 33 tt [1 — 2/ (3 g^) + 5/(33 g'^)] (1 — l/g^)^^; its dependence with g is displayed on Fig. [|. 
We can see that A increases monotonously as g decreases, and that the change is rather steep only for g < 5 (hatched 
region). Otherwise, A stays very close to its asymptotic limit 33 7r « 100. 

Next, we have to relate g with the density n = N/V. We do so by using an intermediate physical parameter, namely 

the screening length k given by = 4 tt — (fi being the chemical potential) . In the case of the Debye-Hiickel 

screening or the Thomas-Fermi model at T = 0, is exactly proportional to n In the latter case, at T 0, 
this remains approximately true provided that n is not too large. The second step is obtained by following Banyai 
and Koch to state that kUx csc 1/g . The value of the proportionality constant depends on how one chooses to 
correlate the Yukawa and Hulthen potentials, but it is always very close to 1 [Q. From these results, we propose the 
approximate relationship 



1 _ n 

g2 ric 



(19) 



where ric is the Mott density, thereby ensuring that no bound state survives at n = ric- For densities n < lO^^nc, 
we get g > 10. Consequently, the influence of screening on the criterion obtained in Ref. |^ is hardly significant. If 
anything, it would lower further the "critical" density above which excitons cannot be considered as perfect bosons. 
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VI. CONCLUSION 



We have shown that screening may be included consistently in calculations pertaining to the departure of excitons 
from purely bosonic behavior. This improvement does not alter substantially a former evaluation, in which screening 
was utterly neglected: Excitons cannot be considered as perfect bosons anymore well below the Mott density (about 
two orders of magnitude). This should not be without consequence on the possible observation of the quite challenging 
Bose condensation of excitons. 
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A more compact expression of Fn is Fn ~ 



where E{r]) = JJ^ (1 + ?7 l^fc 1^) 
v=Q k 
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